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ABSTRACT 


The  nonlinear  potential  equation  that  determines  the  near  aero- 
dynamic and  far  acoustic  fields  of  a hovering  transonic  helicopter 
blade  is  derived.  Properties  of  this  equation  provide  a qualitative 
description  of  characteristic  surfaces  associated  with  the  nonlinear 
flow  field. 

A simple  linearized  acoustic  formula  is  developed;  a nonlineari- 
zation scheme  is  then  applied  to  this  formula,  and  the  linear  and 
nonlinear  acoustic  pressure  profiles,  including  propagating  acoustic 
shock  waves,  are  numerically  determined  and  compared  with  experi- 
mental data. 

A tentative  conclusion  is  that  the  radiated  sound  field  should 
become  measurably  different  when  the  blade  tip  first  punctures  the 
sonic  surface  associated  with  the  rotor  differential  equation.  This 
conjecture  is  suggested  by  the  mathematical  properties  of  the  rotor 
equation  and  recent  experimental  work. 


1.  INTRODUCTION 


It  is  now  fairly  well  established  that  nonlinear  effects  play  an  im- 
portant role  in  the  propagation  of  sound  generated  by  high  speed  helicopter 
blades.  Calculations  of  rotor  blade  acoustic  fields  based  on  linear 
acoustic  theory  have  been  presented  by  several  authors.  Examples  are 
the  work  of  Farassat  (ref.  1),  Hawkings  and  Lowson  (ref.  2),  and  Schmitz 
and  Yu  (ref.  3).  These  calculations  are  in  good  enough  agreement  to  sup- 
port the  conclusion  that  linear  theory  does  not  adequately  describe  recent 
experimental  data.  In  - flight  experimental  data  of  Vause,  Schmitz,  and 
Boxwell  (ref.  4),  and  data  taken  in  an  anechoic  facility  designed  for  hover- 
ing rotor  accoustic  experiments  (ref.  5)  show  wave  form  amplitudes 
that  are  50  to  7 5%  greater  than  amplitudes  based  on  linear  theory.  The 
experimental  data  also  show  a distortion  in  wave  form  that  cannot  be 
accounted  for  by  linear  acoustics.  An  interesting  feature  of  the  data  is 
that  they  show  a wave  profile  that  carries  a pronounced  acoustic  shock 
wave  into  the  far  field.  This  far  field  discontinuity  appears  in  the  simplest 
experimental  arrangement.  For  example,  in  reference  5,  the  model  is  a 
two-bladed  rotor  with  rectangular  planform  , untapered  in  thickness, 
with  aspect  ratio  13.71  andNACA-0012  sections.  The  blade  rotates  in 
hover,  and  some  runs  have  a thrust  coefficient  that  is  just  large  enough 
to  keep  the  blade  wakes  out  of  the  plane  of  rotation.  The  idea  is  to  have 
the  blade  at  near  zero  lift  and  then  effectively  to  measure  an  in-plane 
acoustic  profile  that  is  caused  almost  entirely  by  blade  thickness.  The 
question  is  whether  a linear  theory  of  thickness  sound  predicts  this 
measured  profile  ; the  answer  is  apparently  that  it  does  not. 
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It  is  tempting  to  reach  two  conclusions  here.  Knowing  that  the 
experimentally  determined  acoustic  profile  carries  a shock  wave,  and 
knowing  the  form  of  the  profile  provided  by  linear  theory,  the  conclusions 
would  be  (1)  that  the  experiments  are  displaying  a wave  form  that  is  much 
like  the  profile  found  in  sonic  boom  work,  and  (2)  the  linear  theory  may 
then  be  adjusted,  using  standard  nonlinearization  techniques,  to  give  the 
propagating  shock  embedded  in  the  distorted  profile  shown  in  the  exper- 
imental data.  Both  conclusions  are  approximately  correct.  Clearly,  the 
data  show  that  disturbances  of  an  aerodynamic  nature  on  and  near  the 
blade  surface,  in  its  tip  region,  persist  in  a region  off  the  blade  tip  in 
such  a way  that  a wave  forms;  this  wave  propagates,  rolls  over,  breaks, 
and  forms  a shock  that  moves  off  into  the  far  field.  It  is  also  clear  that 
in  this  picture  there  is  an  important  distinction  between  flow  events  and 
pressure  disturbances  on  and  near  the  blade  tip,  and  pressure  disturb- 
ances at  some  distance  (some  few  chords)  from  the  blade  tip.  .For  contrast, 
consider  a fixed  wing  of  moderate  aspect  ratio,  say  6,  with  rectangular 
planform  and  thickness  ratio  of  about  0.  1,  that  moves  at  constant  velocity 
with  zero  lift  at  Mach  number  of  around  0.  9 or  0.  95.  It  is  well  known 
that  a region  of  embedded  supersonic  flow  occurs  next  to  a portion  of  the 
wing  surface.  Compression  and  expansion  waves  appear  in  this  super- 
sonic region  that  may  be  terminated  by  a shock  wave.  However,  the  wave 
system  is  trapped  between  the  wing  surface,  the  terminating  shock  wave, 
and  a bubble-like  sonic  surface  that  forms  over  a part  of  the  wing.  No 
waves  escape  to  infinity  - there  is  no  boom  or  bang.  The  conterpart  of 
the  fixed  wing  here  is  the  rotor  blade  with  larger  aspect  ratio  - between 
10  and  15  say,  rectangular  planform,  thickness  ratio  around  0.  1,  and 
blade  tip  Mach  number  of  0.  9 - 0.  95.  The  blade  may  be  in  hover  at  zero 
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or  near  zero  lift.  Then  the  flow  field  near  the  blade  surface,  near  the  tip, 
is  very  similar  to  the  field  in  the  supersonic  region  of  the  fixed  wing  . 

The  wing  has  constant  free  stream  velocity,  whereas  the  free  stream 
velocity  of  the  blade  varies  linearly  along  the  span,  increasing  toward 
the  tip.  But  otherwise  a plot  of  pressure  distribution  on  the  rotor  blade 
surface  is  quite  similar  to  the  distribution  of  the  wing,  and  the  location 
and  strength  of  shocks  on  the  two  surfaces  are  also  quite  similar.  The 
two  surfaces  show  similar  aerodynamic  flow  fields  in  their  embedded 
supersonic  regions.  It  is  for  this  reason  that  it  is  convenient  to  classify 
the  fluid  dynamics  near  the  high  speed  rotor  tip  as  nonlinear  and  aero- 
dynamic in  nature,  as  opposed  to  acoustic. 

The  similarities  described  above  appear  in  experimental  data,  and 
also  in  numerical  calculations  based  on  the  small  disturbance  transonic 
potential  flow  equation.  Examples  of  transonic  near  field  calculations 
appear  in  two  papers  by  Caradonna  and  Isom  (refs.  6,  7). 

Any  similarity  between  the  fixed  wing  and  rotor  blade  problems 
disappears  when  the  two  far  fields  are  examined.  Pressure  disturbances 
of  the  transonic  (M.  < 1)  fixed  wing  show  a relatively  rapid  aerodynamic, 
nonacoustic  decay,  as  they  must  if  no  sound  is  produced,  while  disturb- 
ances from  the  rotor  blade  show  a slower,  acoustic  decay  (the  linear 
theory  gives  the  usual  R * geometric  attenuation  of  a wave  profile;  non- 
linear theory  gives  a somewhat  faster  decay).  The  distribution  in  the  far 
fields  is  taken  up  in  Section  2,  where  the  complete  small  disturbance 
differential  equation  for  a hovering  rotor  blade  is  derived.  This  rotor 
differential  equation  neatly  summarizes  the  kinds  of  nonlinearities  that 
affect  the  flow.  The  equation  displays  one  nonlinearity  that  is  effective 
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near  the  blade  surface.  Its  only  role  is  in  the  near  field  aerodynamics. 
One  can  see  by  inspection  of  the  equation  that  this  nonlinearity  becomes 
weak  and  relatively  unimportant  far  from  the  blade.  In  addition,  the 
rotor  equation  shows  a second  nonlinearity  that  does  not  appear  in  sub- 
sonic fixed  wing  work.  This  second  nonlinearity  is  not  effective  on  the 
blade  surface  - it  seems  to  have  nothing  to  do  with  aerodynamics;  its  near 
field  ineffectiveness  is  easily  shown  by  scaling  arguments.  This  non- 
linearity decays  with  distance  from  the  blade  more  slowly  than  does  the 
aerodynamic  one.  It  is  finally  the  effect  that  maintains  distortion  and  a 
shock  wave  in  the  pressure  profile.  It  is  the  acoustic  nonlinearity.  How- 
ever, a mere  inspection  of  the  differential  equation  does  not  show  the 
extent  of  overlap  - in  some  region  off  but  near  the  blade  tip  - of  the  two 
nonlinear  effects  . Experimental  data  (ref.  8)  suggest  that  a shock  provided 
by  aerodynamic  effects  on  the  blade  surface  merges  . somewhere  off  the 
tip,  with  the  more  distinctively  acoustic  shock,  the  latter  discontinuity 
then  propagating  off  into  the  far  field.  Apparently  we  do  not  have  a simple 
boom  problem.  It  is  probably  not  sufficient  first  to  determine  the  linear 
acoustics  of  the  rotor,  and  then  to  adjust  the  linear  result  with  the  usual 
nonlinearization  schemes.  Near  field  effects  seem  to  spoil  this  simple 
program.  Nevertheless,  it  is  this  kind  of  calculation  that  is  taken  up  in 
Sections  3 and  4.  The  recommended  alternative  approach  would  be  a 
full  scale  numerical  attack  on  the  nonlinear  differential  equation  itself. 
Results  from  the  latter  kind  of  calculation  are  not  yet  available. 

We  first  put  linear  theory  in  a very  compact  form  using  transonic 
scaling  arguments  and  approximations.  There  is  nothing  arbitrary  about 
these  initial  approximations  - they  could  be  improved  if  improvement  were 
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required.  The  results  of  Section  3 are  used  to  find  and  plot  the  linear 

wave  profile.  We  then  have  an  F-function  on  which  to  base  a nonlinear 

acoustic  calculation  which  is  taken  up  in  Section  4.  The  relevant  methods 

here  are  all  contained  in  Whitham's  book  (ref.  9).  The  linear  calculation 

gives  an  amplitude  that  is  too  small  in  comparison  with  experimental 

data,  and  this  profile  cannot  distort  as  it  propagates.  The  calculation  of 

Section  4 gives  a distortion  that  has  the  correct  qualitative  appearance, 

and  the  nonlinearization  provides  the  acoustic  shock  wave  in  the  advancing 

profile.  The  adjusted  linear  theory  cannot  give  a shock  that  begins  on  the 

blade  surface,  as  experiment  requires.  At  the  end  of  Section  4 we  modify 

the  nonlinearization  scheme  in  a way  that  is  not  entirely  justified.  It  is 

noted  that  linear  theory  gives  disturbance  pressure  proportional  to 

thickness  ratio  6 to  the  first  power,  whereas  transonic  small  disturbance 

theory  gives  disturbance  pressure  proportional  to  thickness  ratio  to  the 

two-thirds  power.  We  then  simply  rescale  the  linear  pressure  profile 

2/3 

to  make  it  proportional  to  6 . This  modification  affects  the  amplitude 

of  the  wave,  the  distance  from  the  blade  tip  at  which  a shock  wave  first 
appears,  and  the  amount  of  distortion  in  the  profile.  Rescaling  gives  an 
amplitude  that  is  larger  than  the  linear  amplitude  by  a factor  6-1^ 

(which  means  larger  by  a factor  of  2 for  a 12%  thick  blade),  and  the  shock 
moves  in  much  closer  to  the  blade.  While  these  trends  put  the  results  in 
better  agreement  with  experimental  data,  the  method  used  to  obtain  them 
is  not  entirely  satisfactory. 

It  may  be  worthwhile  to  comment  on  the  methods  used  in  this  paper. 
There  seem  to  be  two  quite  different  ways  to  approach  the  high  speed 
helicopter  acoustic  problem.  The  first  method  would  be  based  on  the 


Ffowcs  Williams  - Hawkings  equation  (ref.  10).  The  main  thrust  of  most 
recent  theoretical  research  in  rotor  acoustics  has  been  based  on  some 
use  or  adaptation  of  this  equation.  In  its  most  general  form,  it  is  a con- 
version of  the  equations  of  conservation  of  mass  and  momentum  of  fluid 
dynamics  to  a nonlinear  scalar  integral  equation,  with  boundary  conditions 
on  the  body  surface  accounted  for.  It  is  a generalization  of  Lighthill's  equa- 
tion, the  equation  that  has  been  so  useful  in  jet  noise  work.  Our  linear 
calculation  of  Section  3 is  based  on  the  linear  term  of  the  Ffowcs  Williams  - 
Hawkings  equation  (as  are  other  linear  calculations).  Our  final  form  of 
this  linear  term  is  unlike  any  other  that  has  been  used  to  study  the  rotor 
far  field.  It  would  be  desirable  to  extract  the  essential  nonlinear  features 
of  the  FW-H  equation  and  use  them  in  this  application.  The  main  problem 
is  that  the  Lighthill  stress  tensor  or  its  derivatives  must  be  integrated 
over  all  space  exterior  to  the  body  surface.  Convergence  of  this  integral 
is  very  slow,  and  it  is  not  yet  clear  how  to  approximate  the  stress  tensor 
in  a rational  way,  or  how  sensitive  acoustic  calculations  are  to  approxima- 
tions of  this  nonlinear  term.  This  is  not  to  say  that  the  FW-H  equation 
is  unsuitable  for  the  rotor  problem.  Indeed,  it  seems  unlikely  that  better 
calculations  can  or  should  be  entirely  separated  from  this  kind  of  method. 
The  point  is  that  the  equation  is  still  poorly  understood  in  applications  to 
nonlinear  acoustic  problems  in  infinite  domains  where  wave  profiles  carry 
a shock  wave.  The  equation  has  certainly  not  been  fully  exploited. 

Alternatively,  the  aerodynamic  and  acoustic  field  of  the  rotor  may 
be  formulated  in  terms  of  the  small  disturbance,  nonlinear  partial  dif- 
ferential equation  for  the  velocity  potential.  This  equation  is  discussed 
in  Section  2.  It  seems  likely  that  this  equation  will  play  an  important  role 
in  the  rotor  acoustics  area.  Much  is  known  about  the  numerical  solution  to 
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to  small  disturbance  transonic  potential  equations  in  the  fixed  wing  area. 
The  rotor  potential  equation  has  a most  unusual  feature  that  distinguishes  it 
from  the  fixed  wing  equation.  There  is  a cylindrical  region  centered  on 
the  axis  of  rotation  within  which  the  flow  is  elliptic  and  where  Dirichlet 
boundaiv  condtions  are  applied  on  the  blade  surface.  The  equation  becomes 
hyperbolic  in  the  region  exterior  to  the  cylinder,  and  there  radiation  con- 
ditions must  be  imposed.  The  region  outside  the  cylinder  has  real  charac- 
teristics, and  it  is  these  characteristics  that  carry  off  disturbances  from 
the  blade  to  the  far  field.  These  characteristics  from  an  envelope  (always, 
for  any  non  zero  blade  tip  Mach  number)  and  then  a shock  wave;  it  is  the 
shock  that  has  been  experimentally  measured.  The  differential  equation 
is  used  in  Section  4 to  find  corrections  to  the  constant  ambient  sound  speed 
in  the  far  field,  and  then  to  improve  the  characteristics  of  the  linear  rotor 
potential  equation.  This  procedure  is  called  nonlinearization.  It  is  now 
well  known  from  sonic  boom  work. 
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2.  THE  EQUATION  OF  MOTION  AND  ITS  NONLINEARITIES 


The  simplifying  assumptions  for  flow  over  a transonic  helicopter 
blade  are  essentially  the  same  as  those  for  transonic  flow  over  a fixed 
wing:  the  gas  has  constant  specific  heats,  any  shock  waves  are  weak  and 
entropy  jumps  across  them  may  be  neglected,  the  disturbance  velocity 
is  irrotational  and  therefore  has  a potential,  and  only  second  order  non- 
linearities  need  be  retained  in  the  potential  equation.  With  these  as- 
sumptions the  differential  equation  of  motion  is 


V0t,  = a 2V  V 


(2.  1) 


The  required  approximation  to  the  adiabatic  sound  speed  squared  is 


a2  = aQ  - (Y-l  )0t,  . (2.  2) 

with  aQ  the  sound  speed  in  undisturbed  air  and  u = V0.  Equations  (2.  1) 
and  (2.  2)  are  correct  for  a coordinate  system  fixed  in  undisturbed  air. 

We  consider  flow  over  a hovering  helicopter  blade;  that  is,  flow  over  a 
blade  that  rotates  at  constant  angular  velocity  w and  that  has  no  forward 
motion  (zero  advance  ratio).  Choose  a system  of  cylindrical  coordinates 
that  rotates  with  the  blade  system  and  relative  to  which  the  blades  are  at 
rest.  Let  these  coordinates  be  r,  6,  z,  with  positive  z pointing  in  the 
direction  of  the  angular  velocity  vector  w and  perpendicular  to  the  plane 
of  rotation.  Figure  1 shows  a planview  of  two  rectangular  blades.  The 
picture  is  simplified  by  ignoring  the  geometry  of  the  blades  near  the 
rotor  hub  — there  is  no  cut-out.  Also  shown  are  Cartesian  coordinates 
x',  y'  fixed  in  undisturbed  air,  and  coordinates  x,  y along  the  blade 
chord  and  blade  span,  and  rotating  with  the  blade.  An  observer  at 
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point  O in  the  plane  of  rotation  has  radial  position  r,  polar  angle  0' 
relative  to  the  nonrotating  y'  axis,  and  angle  0 relative  to  the  rotating 
span  axis  y.  The  cylindrical  coordinate  systems  are  related  by 


r = r' , 

0 = 0'  -ut'  , 

Z - z'  , 


and  the  two  time  scales  are  the  same,  t'  = t.  The  primed  system  is  the 
inertial  one,  and  the  unprimed  system  is  blade-fixed.  The  potential  <£ 
for  a hovering  rotor  is  a function  only  of  the  variables,  r,  0*-«t,  and  z. 
Time  derivatives  in  (2.  1 ) transform  to 

2 

<pt,  = x r • r <p  = 4>t,t,  = w (pee . 


Using  these  results  and  some  rearrangement  of  terms,  equations  (2. 
and  (2.  2)  combine  to  give 


a 

-f-(Y 

r 


1 * ■% 
r 


4>ee  - ^Ae  - 


l) 


a^  + (Y-D 


(2.3) 


In  all  that  follows,  we  need  to  consider  only  the  flow  over  one  blade 
of  the  two-bladed  system.  The  reason  is  that  a single  blade  has  large 
aspect  ratio,  generally  greater  than  10,  and  the  main  concern  is  with 
flow  events  near  a blade  tip  where  the  flow  may  be  transonic.  A good 
approximation  is  that  flows  near  the  two  blade  tips  are  essentially  in- 
dependent of  each  other.  When  we  come  to  acoustic  considerations,  the 
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most  intense  part  of  the  sound  field  is  generated  by  a single  tip,  the  one 
that  is  rotating  toward  an  observer.  We  choose  the  blade  on  which 
y > 0. 

Boundary  conditions  for  (2.  3)  are  required.  The  main  problems 
of  rotor  aerodynamics  and  acoustics  are  encountered  with  the  simplest 
blade  geometry:  rectangular  planform,  no  taper  in  thickness  along  the 
span,  and  constant  section  geometry.  Figure  2 shows  this  simple 
geometry  in  blade  fixed  coordinates.  A blade  has  chord  i and  span  L. 
We  take  the  y-axis  as  the  blade  span  axis  running  through  section  quar- 
ter chord  points.  The  leading  edge  of  the  y > 0 blade  is  x = -i/4  and  the 
trailing  edge  is  at  x = 3i/4.  This  disposition  of  coordinate  axes  is  tem- 
porarily convenient.  It  will  be  more  convenient  later  on  to  place  the 
y-axis  along  the  leading  edge  of  the  blade.  We  take  blade  section  geo- 
metry to  be 

G(x,  y,  z)  = z - 6 if^|  + -0  = 0,  (0  < y < L) 

where  6 is  the  thickness  ratio.  Now  let  V be  the  velocity  of  air  rela- 
tive to  the  blade,  so  that 

V = - u x r + V0. 

The  condition  of  tangency  at  the  blade  surface  is  V • VG  = 0 on  G = 0. 

A linearized  boundary  condition  with  the  mean  surface  approximation  is 

J£=«y«f  on  z a 0,  (2.4) 

the  dot  over  f indicating  a derivative  with  respect  to  x/jt. 

The  behavior  of  a solution  to  equations  {2.  3)  and  (2.  4)  in  the  far 
field  must  be  specified.  The  required  behavior  results  from  imposing  a 
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radiation  condition  for  r -*  oo  . The  form  that  this  condition  takes  is 
most  easily  found  by  first  transforming  (2.  3)  to  spherical  coordinates 
R,  a , 9.  Let  R be  the  position  vector  from  the  center  of  rotation  to 
a variable  point,  with  |r|  = R,  a the  angle  that  R makes  with  the  posi- 
tive z-axis,  and  9 the  same  in -plane  polar  angle  as  in  Figure  1.  Equa- 
tion (2.  3)  in  this  coordinate  system  is 
2 . 


2 

w 


u)( 


2 5 " (Y  "*■  1 ) 2 2 

R sin  or  R sin  & 


%9  ' 2w<^R^R 9 “ 2 r2  Wae 


ao  “ 


^RR  + R ^R  + „2  . da  (<#,a8inQ') 
R sines 


. (2.5) 


It  is  of  interest  to  examine  the  form  this  equation  takes  for  a linearized 
problem.  Neglecting  nonlinear  terms  in  (2.  5)  gives  the  simple  equation 
2 


2 °o  \ , 2 

"2.2  )(t>99  ~ ao 

R sin  a 


^RR  + R ^R  + 2 . aa  (<Pa9lna) 

K sin  a 


(2.6) 


Equation  (2.  6)  has  the  unusual  property  that  it  is  elliptic  for 
2 

2 . a- 

w < 


D 2 . 2 

R sin  or 


and  hyperbolic  for 
2 

2 . ao 

“ > _2  72  ' 

R sin  a 

In  cylindrical  coordinates,  linearization  of  equations  (2.  3)  gives  elliptic 
behavior  for 

2 *o 

« < -J. 

r 


« ' .ju  t 
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and  hyperbolic  behavior  for 
2 

2 ^ ao 

« > 2 ' 
r 

The  conclusion  is  that  in  linear  theory  there  is  a cylinder  of  radius 
r = aQ/«  inside  which  the  flow  is  elliptic  and  outside  which  the  flow 
is  hyperbolic.  We  also  conclude  that  no  waves  or  wave  propagation 
may  appear  inside  this  cylinder,  while  wave  systems  and  wave  propa- 
gation may  occur  outside  it.  The  cylinder  r = aQ/<«o  is  called  the  sonic 
cylinder. 

The  tip  Mach  number  M of  a blade  is  defined  as  M = «L/aQ, 
where  L is  the  blade  radius.  The  sonic  cylinder  is  then  defined  by 

r _1_ 

L “ M- 

The  above  comments  refer  to  a linearized  description  of  the  flow,  but 

it  is  evident  that  the  main  difficulties  in  calculating  flow  over  a transonic 

rotor  are  intimately  related  to  the  fact  that  for  M near  one,  the  sonic 

cylinder  lies  very  close  to  the  blade  tip.  Referring  to  the  nonlinear 

equation  (2.  3),  we  will  have  elliptic  and  wave -free  behavior  for 
a2 

j - (y  + 1 ) -Tj-  <()q  < 0,  (2.  7a) 


and  hyperbolic  and  wave -like  behavior  for 
2 


W*  - “ (Y  + 1 ) <f>Q  ^ 


0. 


(2.7b) 


t 

Of  course  0^  is  not  known  a priori,  and  so  the  shape  of  the  surface  de- 
fined by  setting  the  left-hand  side  of  (2.  7)  equal  to  zero  is  not  known 
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until  a solution  is  found.  Nevertheless,  we  may  conclude  from  (2.  7) 
that  the  sonic  cylinder  of  linearized  theory  is  replaced  in  a nonlinear 
theory  by  a kind  of  deformed  cylindrical  surface.  If,  for  example,  r 
and  6 are  held  fixed  while  z - + oo  in  (2.  7),  then  - 0 asymptotically, 
and  the  surface  that  separates  elliptic  from  hyperbolic  behavior  returns 
(asymptotically)  to  the  cylinder  r = aQ/w . One  may  then  say  that  the 
sonic  surface  defined  by  (2.  7)  will  depart  from  true  cylindrical  form 
m a region  surrounding  the  blade  tip,  where  the  surface  may  have  bulges 
toward  or  away  from  the  blade  tip.  The  surface  defined  by  setting  the 
left-hand  sides  of  (2.  7)  equal  to  zero  is  called  the  sonic  surface. 

It  is  tempting  to  make  the  following  conjecture:  Peak  or  high 
intensity  sound  generated  by  a transonic  helicoptor  blade  first  occurs 
when  the  blade  tip  penetrates  the  sonic  surface.  The  idea  is  that  the 
differential  equation  (2.  3)  has  a region  of  hyperbolic  behavior  in  which 
there  are  real  characteristics  when  the  inequality  (2.  7b)  holds.  Such 
a region  exists  for  an^  non  zero  tip  Mach  number,  no  matter  how  small 
or  large  M may  be.  If  M is  very  small,  the  space  in  which  character- 
istics exist  is  at  some  great  distance  (roughly  1 /M)  from  the  center  of 
rotation.  The  characteristics  do  not  touch  the  blade  surface.  As  M ♦ 1 
from  below,  the  sonic  surface  draws  nearer  to  the  blade  tip.  The  far- 
ther this  surface  is  from  the  tip,  the  more  the  shape  is  like  a circular 
cylinder.  The  nearer  the  surface  is  to  the  tip,  the  more  it  will  be  de- 
formed near  the  tip.  There  must  then  be  some  tip  Mach  number  at 
which  the  sonic  surface  first  comes  in  contact  with  the  blade  tip.  If  M 
is  increased  beyond  this  lower  critical  Mach  number,  a portion  of  the 
blade  tip  will  protrude  through  the  sonic  surface.  Real  characteristics 
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on  which  disturbances  propagate  are  then  in  contact  with  the  blade  sur- 
face, and  these  characteristics  escape  to  infinity  (i.  e. , to  the  far  field). 

There  is  as  yet  neither  computational  nor  experimental  evidence 
to  support  or  disprove  these  qualitative  remarks.  Recent  experimental 
work  by  Schmitz  et  al.  (ref.  3 ) and  the  recent  report  by  Yu  et  al.  (ref. 

8 ) suggest  that  our  guess  (conjecture)  and  the  above  remarks  may  be 
interesting. 

We  return  to  equation  (2.6),  the  linearized  potential  equation  in 

spherical  coordinates.  Let  the  meridian  angle  a be  fixed  and  not  zero, 

2 

let  R — oo  , and  neglect  terms  divided  by  R . Then  (2.6)  becomes 


which  is  also 

«2(R0)0e  = a2(R0)RR.  (2.  8) 

This  result  shows  simple  behavior  of  solutions  to  (2.6)  in  the  far  field, 
for  (2.  8)  is  just  a wave  equation  with  solutions 

R(f)  = Fj  (uR  + aQ0  ) + FglwR  - a q0  ).  (2.  9) 

The  characteristic  surfaces  of  (2.  8)  are 

uR  + a 0 = c,,  wR  - a 0 = c„ 

or  o 2 

where  c^  and  Cg  are  constants.  The  trace  of  these  surfaces  in  the 
plane  of  rotation  z = 0 are  the  Archimedean  spirals 

wr  - a 0 = c . 

O 4 


ur  + a 0 = c, , 
o r 
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Only  one  of  the  functions  and  Fg  in  (2.  9)  represents  the  solution  0 
for  large  R.  The  choice  is  based  on  the  fact  that  small  disturbances 
propagate  in  a "downstream''  direction.  For  an  observer  fixed  relative 
to  the  blade,  the  free  streamlines  are  circles  with  centers  on  the  axis 
of  rotation.  If  the  blade  rotates  in  a counterclockwise  direction  relative 
to  undisturbed  air,  then  an  observer  fixed  on  the  blade  would  see  an  un- 
disturbed air  motion  in  a clockwise  direction,  as  in  Figure  3.  The 
figure  shows  arcs  AB  and  AC  of  the  traces  in  the  plane  z = 0 of  the  far 
field  characteristics  c^  and  Cg.  The  c^  characteristic  is  the  one  that 
carries  a disturbance  at  large  R in  the  downstream  direction  from  A to 
B.  Hence  Fj  is  the  proper  choice  and  we  have 

R0~F(wR  + aQ0,  a-)  (R -*•(»,  0 < a < x/2).  (2.10) 

A subscript  on  F is  no  longer  needed. 

The  linearized  approximation  (2.  10)  to  the  acoustic  potential  shows 
that  wave  propagation,  in  this  approximation,  is  without  distortion  and 
subject  only  to  geometric  attenuation.  This  behavior  is  a result  of  as- 
suming a constant  speed  of  sound  in  the  linear  theory.  A more  accurate 
description  of  far  field  wave  motion  can  be  based  on  the  linear  result 
(2.  10)  with  a nonlinear  correction  to  the  local  sound  speed.  This  tech- 
nique results  in  a slightly  faster  than  R * attenuation,  produces  dis- 
tortion in  the  wave  form  F,  and  finally  results  in  the  formation  of  an 
acoustic  shock  wave.  This  method  is  applied  and  discussed  in  Section  4. 

It  is  a fairly  straightforward  matter  to  nonlinearize  the  approxi- 
mation (2.  10)  by  correcting  the  far  field  sound  speed.  The  main  prob- 
lem here  is  to  find  the  right  wave  form  F to  nonlineariae.  An  F based 
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on  the  strictly  linear  theory  is  one  candidate,  but  use  of  a linear  F im- 
plies that  the  flow  field  near  the  blade  surface  is,  in  some  sense,  well 
approximated  by  a linear  theory,  and  hence  that  nonlinear  corrections 
are  required  only  in  the  far  field.  This  kind  of  problem  — linear  near 
field,  nonlinear  far  field  — often  appears  in  the  simplest  sonic  boom 
work.  An  unusual  feature  of  the  hovering  rotor  differential  equation 
(2.  3)  is  that  it  has  two  kinds  of  nonlinearities.  The  first  one  appears 
in  a region  surrounding  the  blade  tip.  It  is  an  aerodynamic  nonlinearity 
and  is  associated  with  the  transonic  nature  of  flow  in  the  blade  tip  region. 
If  an  observer  point  is  near  or  in  the  plane  of  rotation,  and  the  radial 
coordinate  r is  nearly  equal  to  the  blade  radius  L,  then  equation  (2.  3) 


is  well  approximated  by 
2 

2 S , . . i 

“ “ “2  " (V  + 1)  ~2 

r r 


= A (<f)  + <f>). 

n “rr 


zz 


(2. 11) 


The  problem  posed  by  (2.  11)  and  the  tangency  condition  is  roughly  simi- 
lar to  the  nonlinear  transonic  problem  for  a fixed  wing.  However,  any 
similarity  between  the  fixed  wing  and  rotor  problems  disappears  in  the 
far  field.  The  far  field  rotor  problem  is  best  described  by  switching 
from  the  cylindrical  equation  (2.  3)  to  the  corresponding  spherical  form 
(2.  5).  In  (2.  5),  let  R -*•  oo  while  the  meridian  angle  a is  fixed  and  not 
zero.  Then  (2.  5)  reduces  to 


2 , 


GG 


- 2“<t>R<t>Re  = 


aQ  + (y  - 1 ) w , 


(2.  12) 


A slow  variation  of  local  sound  speed  for  large  R is  embedded  in  this 
equation.  The  quadratic  terms  here  are  properly  called  acoustic  non- 
linearities.  They  are  small  perturbations  to  an  otherwise  linear 
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problem,  as  opposed  to  the  nonlinearity  in  equation  (2.11)  which  is  not 
small  relative  to  other  terms  in  that  equation. 

This  schematic  division  of  the  flow  field  into  an  inner  aerodynamic 
tip  solution  and  an  outer  acoustic  field  is  not  essential.  In  principle,  a 
numerical  solution  could  be  based  on  the  complete  differential  equation 
(2.  3),  with  a surface  tangency  condition  and  an  outer  numerical  boundary 
condition  that  simulates  a radiation  condition.  Acoustic  results  from 
this  kind  of  calculation  are  not  yet  available,  although  the  near  field 
aerodynamic  problem  has  been  explored  by  Caradonna  and  Isom  (refs. 

6 and  7 ). 

The  alternative  course  is  to  base  the  outer,  nonlinear  acoustic 
solution  on  the  linear  solution  (2.  10)  and  a corrected  sound  speed  that 
is  extracted  from  (2.  12).  This  is  essentially  a sonic  boom  calculation. 
The  simplest  linear  problem  is  the  one  that  accounts  for  sound  due  to 
blade  thickness  alone.  An  explicit  solution  is  available.  In  the  follow- 
ing section  we  shall  put  that  solution  in  its  limiting  transonic  form.  The 
result  is  a very  compact  formula  that  is  suitable  for  a nonlinearization 
scheme  that  is  required  in  an  acoustic  shock  calculation. 
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3.  LINEAR  THEORY  IN  THE  TRANSONIC  LIMIT 


The  disturbance  pressure  p in  the  linear  theory  of  potential  flow 
over  a rigid  body  is  given  by  the  well  known  formula 


i 9 f n dS 

4,rp  ‘ po  at  J i -mr  R ' 


(3.  1) 


The  symbols  have  their  usual  meaning.  On  the  right-hand  side,  pQ  is 
constant  ambient  density.  If  v is  the  velocity  of  a point  on  the  body 
surface  relative  to  stationary  air,  then  v^  is  the  projection  of  v along 
the  outward  normal  to  the  body  surface;  MR  is  the  radiation  Mach  num- 
ber, and  the  script  or  fancy  JZ  is 

/?=|R-R0I,  (3.2) 


where  R is  the  observer's  position  vector  and  R is  the  position  vector 

o 

to  an  integration  point  on  the  body  surface.  The  integral  is  over  the 
body  surface  and  all  terms  in  the  integrand  are  evaluated  at  retarded 
time  t . This  time  is  defined  by  the  retarded  time  equation 

ao(t-r)=/?,  (3.3) 


where  a is  the  constant  ambient  adiabatic  sound  speed, 
o 

For  a hovering  helicopter  blade,  it  is  convenient  to  consider  the 
problem  of  determining  p as  defined  in  blade  fixed  coordinates.  An 
observer  fixed  in  undisturbed  air  would  then  appear  to  rotate  relative 
to  the  blade  in  a clockwise  motion  about  the  blade  center,  if  the  blade 
motion  is  counterclockwise  relative  to  stationary  air.  If  we  let 
n ~ (rjj,  r)2*  13)  be  Cartesian  blade  fixed  coordinates,  then  Rq  * R^q)  in 
(3.  2)  and  the  observer  position  vector  in  (3.  2)  is  a function  of  retarded 
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time  t,  so  we  have 

/?=  | R(t ) - R0(n)l. 

A 

The  unit  radiation  vector  R is  defined  as 


R = 


|R(r)  - Ro(n)| 


(3.4) 


• — ► A 

In  the  far  field,  as  | R — » , this  vector  is  denoted  by  r instead  of  Rs 


^ _ R(t) 

r - . 


(3.5) 


Let  i,  j,  k be  unit  vectors  along  the  rjj,  rjg'  ^3  axes.  Let  9'  be  a cylin- 
drical polar  angle  in  the  plane  of  rotation,  measured  positive  counter- 
clockwise from  a fixed  axis,  and  let  \p  be  the  angle  that  the  ob- 
server  position  vector  R makes  with  the  plane  of  rotation  so  that 
R • k > 0 for  0 < (/-  < -n  /2.  Then  r is 

A A A A 

r = - i cos  ip  sin(0'  -<jt)  + j cos  \p  cos  (O'  -ut)  + k sin  \p  . (3.6) 

Notation  is  mixed  here.  The  unit  coordinate  vectors  in  (3.6)  are  fixed 
relative  to  the  blade,  while  6 1 is  the  observer  in-plane  angle  measured 
from  a reference  axis  fixed  in  undisturbed  air.  The  quantity  0'  -wt  is 
therefore  the  observer's  retarded  angle  relative  to  the  blade  span  axis 
(relative  to  the  blade,  the  observer  does  the  rotating). 

In  the  far  field,  formula  (3.  1 ) for  p becomes  the  formula  for 
linearized  acoustic  pressure. 

The  Doppler  factor  (l  - M^)’*  appears  so  often  that  we  will  denote  it 
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by  a special  symbol  D; 


D s 


1 


1 -Mr’ 


so  that  acoustic  pressure  p is 

4?rRP  = Po!  /vnDdS< 


(3.  7) 


(3.  8) 


Now  it  is  fairly  easy  to  prove  that  the  acoustic  pressure  p is 
4*RP  = PQao  -j|  / r • n D dS,  (3.9) 

A 

where  n is  the  blade  outward  pointing  unit  normal  vector.  Formulas 
(3.  8)  and  (3.  9)  are  acoustically  or  radiatively  identical.  The  proof  of 
this  equivalence  is  given  by  Isom  in  reference  (II).  The  formula  (3.  9) 
is  more  convenient  for  our  purposes  than  its  equivalent  and  better 
known  form  (3.  8). 

We  now  write  out  the  integral  (3.  9)  in  more  detail.  The  body  is  a 
blade  with  rectangular  planform  and  constant  section  geometry,  un- 
tapered in  thickness.  Let  S.  be  the  blade  chord  and  L its  span,  the 
span  being  the  distance  from  the  center  of  rotation  to  the  blade  tip. 
Coordinate  rj  j is  along  the  chord,  rjg  along  the  span  axis,  and  is 
normal  to  the  blade  mean  surface  and  points  in  the  direction  of  w. 

Section  geometry  is  defined  by 

J?3  = + b). 

The  nondimens ional  number  b is  the  fraction  of  the  chord  from  the  blade 
leading  edge  to  the  blade  span  axis.  This  axis  could  be  taken  through 
the  section  quarter-chord  points;  however,  its  precise  location  does  not 
play  an  important  role  in  acoustic  calculations;  6 is  the  section  thickness 
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ratio.  Define  nondimens ional  blade  coordinates  by 


. r 1?  r Hi 

*1  ~ l ’ ' L ' g3 " I 


Section  geometry  in  these  coordinates  is 


? 3 = 6f(?j  + b). 

A 

The  unit  vector  n normal  to  the  blade  surface  is  then 

— * * ■ o . 9 ■!  /2 

n = (+  k - i 6 f )(  1 + 6^  f*) 


The  plus  sign  on  k denotes  the  blade  upper  surface  0 and  the  minus 

sign  the  lower  surface  ^ 0*  The  projection  of  the  area  element  dS 
onto  the  blade  mean  surface  is  dSQ,  with 


dS  = (1  + 


2 f 2 y2  dS 
6 f ) c 


and 


dSQ=  dni  dT,2=  iLdJjdCj. 

With  these  definitions,  the  formula  (3.  9)  for  acoustic  pressure  p is 

4»Rps  2pQaQf  L ^6  f cos^  sin(0'-tjT)  + sin^/ |d  d?  j d^g*  (3.10) 

A mean  surface  approximation  has  not  yet  been  applied  to  the  integrand 
of  (3.  10);  the  integrand  in  (3.  10)  is  evaluated  on  the  exact  surface  distri- 
bution and  not  on  = 0.  We  are,  however,  assuming  a symmetrical 
section  geometry,  and  the  factor  of  two  in  (3.  10)  accounts  for  the  upper 
and  lower  surfaces.  Thus  the  integral  need  be  taken  only  over  the  upper 


surface, 
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When  the  mean  surface  approximation  is  applied  to  (3.  10),  it  is 
necessary  to  expand  t and  the  Doppler  factor  D as  functions  of  the 
normal  coordinate  $3.  Both  r and  D generally  depend  on  ^ and 

? y but  we  are  here  concerned  only  with  the  dependence,  so  the  re- 
quired expansions  take  the  form 

t<53>=  t(0)  + ^ 6f +...  (3.Ha) 

D(?3)=  D(0)  + -^i6f  +...  (3.11b) 


Function  f appears  on  the  right  because  ?3  is  evaluated  at  ?3  = 6 f , i.  e. , 
on  the  blade  surface.  It  is  tedious  to  apply  the  expansions  (3.  11)  to 
the  integrand  of  (3.  10),  but  the  process  is  straightforward.  The  result 
is 


4*Rp  = 


2poVL6  a 


cos  \p 


J f sin(0'  -wt)D  dl^dSg, 


(3.12) 


where  now  and  henceforth  it  is  understood  that  the  integrand  is  evaluated 
on  ?3  = 0.  The  double  integration  is  from  blade  leading  edge  to  trailing 
edge  and  blade  root  ?2  = 0 to  blade  tip  $2  = 1. 

We  next  introduce  a transformation  in  (3.  12)  that  will  permit  one 
integration  by  parts  in  the  spanwise  direction,  thereby  reducing  the 
formula  p to  a single  integral.  The  result  will  be  that  p is  expressed 
as  a single  integral  over  the  blade  tip  ?2  = *•  The  transformation  pro- 
ceeds as  follows.  It  is  well  known  that  the  gradient  of  retarded  time  r 
with  respect  to  dimensional  coordinates  rjj.  r)2*  ^3  ls 

A 

y ^ t » JL  d . 
n o 
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Here  R is  the  unit  radiation  vector  defined  in  (3.  4).  In  the  far  field 
R — » We  write 


V.t=-D, 
n o 


with  r given  by  (3.  5).  The  derivative  of  r along  the  blade  span  axis 


?2  iS 


9r  A a a Q jn 

= j • V_T=j  • r — = cos  0 cos(0'  -wr)  — . 


dr) 


We  solve  this  equation  for  the  Doppler  factor  D to  get 
a 


D = 


1 


9t 


cos  ip  cos(0'-ut)  dr). 


The  formula  (3.  12)  for  p becomes 

p a2f  L6 
r /\  « 


4irRp  = ° f ’i  tan(0'  -ut)  d?  d£ 


cos  ip 
Substitute  the  identity 


dr). 


(3.  12a) 


“a~log  |cos(0' -wr)|  = tan(9' -ut)^- 
w vrjg  ot]2 


into  (3.  1 2a)  to  get 

p a li  . _ 

4»Rp=  ? A Ji  log|cos(0'-(jT)|dr  dCr 

u cos  0 s2  i ^ 


The  section  slope  distribution  f depends  only  on  so  the  integra- 

1 a 

tion  may  be  carried  out  to  give 


p a 2t6  _ 

4*Rp  = -°-^2  /f  log 

fa)  COS 


cos(9'-ti)T+) 

COS(0'-<JT  ) ^1* 


(3.13) 


Notation  t+  means  evaluate  retarded  time  t at  chordwise  points  on  the 
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blade  tip  section  = 1 > T means  evaluate  t on  chordwise  points  on  the 
blade  root  section  ?2  = °*  The  integral  in  (3.1  3)  is  from  leading  edge  to 
trailing  edge,  on  the  tip  and  root  sections. 

Effective  use  of  equation  (3.  13)  requires  some  tailoring  of  the  re- 
tarded time  equation  (3.  3).  In  what  follows,  approximations  are  based 
on  the  assumption  that  the  blade  aspect  ratio  is  large  and  the  tip  Mach 
number  is  near  one.  For  a single  rectangular  blade  with  chord  4 and 
span  L,  aspect  ratio  is  L/4  . We  define  the  small  parameter  « as  the 

e = 4/L  . 

Approximations  in  the  retarded  time  equation  will  be  based  on  the  as- 
sumption that 


This  approximation  is  analogous  to  the  usual  transonic  fixed  wing  approxi- 
mation that  thickness  ratio  to  the  two -thirds  power  is  small.  The  analogy 

results  because  here  a transonic  acoustic  similarity  parameter  will  ap- 

2/3  2 /3 

pear  that  is  based  one  ' rather  than  6 

We  first  note  that  peak  or  high  intensity  sound  occurs  in  the  far 
field  when  the  quantity  0-ur  in  equation  (3.  13)  is  near  rr/2.  This  result 
corresponds  to  the  fact  that  an  observer  will  measure  peak  sound  when 
his  line  of  sight  is  nearly  perpendicular  to  the  blade  span  axis,  and  the 
blade  is  rotating  toward  him.  We  therefore  define  a retarded  angle  f 
by 

0'  -wt  = ir/2  - € , 


where  5 will  be  small  in  a sense  to  be  determined.  The  formula  (3.  13) 
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for  p in  terms  of  ? is 

.2,, 


4jrRp 


2, 

(J  COS  if/ 


at 


fi  log 

sin  £ + 

sin  £ 

” 

dS 


i * 


(3.14) 


The  plus  and  minus  signs  on  retarded  angle  f have  the  same  signifi 
cance  as  the  subscripted  signs  on  retarded  time  r. 

The  retarded  time  equation  (3.  3)  may  be  written 

9 ^ _ 2 i/2 

ut  -UT--  -^(R2  - 2R-  Ro  + Rc)  • < 

o 

With  R restricted  to  points  on  the  blade  surface  and  R large,  the 

o 

usual  acoustic  approximation  to  (3.  15)  is 


ut 


uR  _u 
ao  'a 


(3.16) 


T,et  R lie  in  the  blade  mean  surface 
o 


Then  (3.  16)  is 


wt  _ ur  = — + — rj.  cos  if/  sin(0'-uT)  - rj2  cos  0 cos(0'-ut). 
ao  ao  1 o 

Nondimens ionalize  „ l and  r?2  by  £ and  L,  use  the  definition  e = £ /L 
and  set  0'-ur  = ir/2  - ? to  get  the  retarded  angle  equation 


£ _ i.  _ - M e ? | COS  if/  cos  ? + M?  2 cos  ip  sin  5 . 

Take  ? to  be  small  and  note  that  the  term  in  cos  f is  already  of  order  e , 
so  that  cos  C may  be  set  equal  to  one  there.  The  term  in  sin  ? requires 
a two-term  Taylor  series  approximation  when  M cos  ip  is  near  one  and 
the  observer  is  in  or  near  the  plane  of  rotation.  We  then  get  an  alge- 
braic equation  for 

5 « | - (e'  -wt  + -Me?!  cos  if/  + M*  2?  cos  ip  - ^ 5 2?  3 cos  if/  . 
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Rear  ranging  terms  gives 

(—■  ?2  cos  + ( 1 - M?2  cos  - (o'-wt  + - M e cos 

(3.17) 

The  retarded  angle  ? must  be  evaluated  on  the  tip  section  where 

£ x l and  5=5.  and  at  the  root  where  ?2  s 0 an<^  5 = ?_•  ^e  tip  and 

2 • 

root  retarded  angle  equations  are  then 

Og-  cos  <//)?++  (1  -Mcos*«+*f -(e'-wt  + ^-McSj  cos  ^ 

(3.18) 

e = ■£  - e'-ut+—  -Mcr  cos  ip.  (3.19) 

- 2 ao  1 

The  cubic  term  in  § on  the  left  of  (3.  18)  is  negligible  when  the  tip  Mach 
number  M is  small,  or  when  the  observer  meridian  angle  ^ is  near 
?r/2,  that  is,  when  the  observer  is  close  to  the  axis  of  rotation.  When 
either  or  both  of  these  conditions  hold,  the  entire  blade  surface  becomes 
acoustically  compact.  On  the  other  hand,  when  M is  near  one  and  the 
observer  is  in  or  near  the  plane  of  rotation  so  that  the  meridian  angle  \p 
is  near  zero,  the  factor  1 - M cos  ^ that  multiplies  ? in  (3.  18)  becomes 
small  and  the  cubic  term  is  no  longer  relatively  small.  In  this  case, 
the  outer  portion  of  the  blade  surface  is  not  acoustically  compact.  The 
presence  of  the  cubic  term  reflects  this  noncompactness  and  the  im- 
portance of  time  retardation.  The  simple  equation  for  results  be- 
cause acoustic  source  points  near  the  center  of  rotation  are  always 
compact;  sound  from  these  low  speed  poinds  has  no  significant  Doppler 
shift. 

The  integral  formula  (3.  14)  requires 
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sin  ? 5 ' 


(3.  20) 


only  the  first  approximations  to  sin  £ for  small  ? being  needed  here.  The 
right-hand  side  of  (3.  20)  is  evaluated  using  (3.  18)  and  (3.  19)  to  give 


1 


§-  cos  0)  ?2  + 1 - M COS  Ip 


(3.21) 


The  final  step  in  these  approximations  is  to  transonically  scale  the  £ + 
retarded  angle  equation  (3.18).  Define  a scaled  characteristic  variable 
fl  by 


tt/2  - (O'  -ut  + {jR/aQ) 


(3.22) 


e M cos  \p  * 

and  a scaled,  retarded  angle  3 by 

e 5 ("-6(1  - M cos  I//)]  1 7/2 

[ M cos  ip  J p ' 

With  these  definitions,  the  ? retarded  angle  equation  (3.  18)  on  the  blade 
tip  becomes 

1 3/2 


3 1 

3 + 3=  j- 

•JH 


M cos  \p 
1 - M cos  \p 


e(n-  ?t). 


(3.23) 


Finally,  we  define  a transonic  acoustic  similarity  parameter  by 


K . -L  f.  M cos  (p  "I 
VT  L1  - M co» 


3/2 


e . 


The  scaled  retarded  angle  equation  (3.  23)  reduces  to 


(3.24) 


3 + 3=  K(n-  k{). 


(3.25) 


X (it 
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If  the  similarity  parameter  K is  small,  then  the  entire  rotor  blade  sur- 
face is  acoustically  compact.  If  K is  0(1)  or  large,  points  on  the  blade 
surface  are  not  compact  and  there  is  a significant  variation  of  retarded 
angle  (and  time)  over  the  surface. 

Now  replace  the  ratio  in  equation  (3.  21 ) by 


? 


_+  1 1 

? 1 - M cos  [p  ^2  + j 


(3.26) 


The  right-hand  side  of  (3.  26)  replaces  sin  £+/sin  ? in  the  integral  for- 
mula (3.  14)  for  p.  Observe  that  when  the  logarithm  of  5 /?  is  taken, 
the  constant  factor  (1  - M cos  ^)~*  in  (3.  26)  may  be  ignored  because  it 
gives  zero  contribution  in  the  integral  (3.  14).  This  is  because  the  lead- 
ing and  trailing  edges  of  the  blade  lie  in  the  plane  of  rotation,  and  the 

i 

thickness  distribution  is  zero  at  these  two  points. 

The  formula  (3.  14)  for  acoustic  pressure  is  now 


4jrRp 


2<V§*S  8 

o 

ucos  </< 


■ft  ft  log(/3Z  + 1)  d^. 


(3.27) 


We  next  replace  the  t-derivative  by  a characteristic  variable  ^-derivative, 
the  5l(t)  relation  being  (3.  22).  Then  (3.  27)  becomes 


4,Rp.  . 2po,^_S_  i fi  logO2  ♦ Ddt, 


(3.28) 


M cos  0 


The  slope  distribution  f , with  its  singular  behavior  at  a round  leading 
edge,  may  be  eliminated  in  favor  of  the  thickness  distribution  f by  an 
integration  by  parts.  This  requires 


* 


K 


3/3Z  + 1 
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from  (3.  25).  Also,  the  O-derivative  in  (3.  28)  may  be  taken  inside  the 
integral.  This  requires 


K 

9 ’ 

3/r  + i 


using  (3.25)  again.  The  last  step  is  to  move  the  origin  of  the  integra- 
tion variable  jjj  to  the  blade  leading  edge.  Let 

X * C,  + b. 


Our  formula  for  acoustic  pressure  p is  finally 


P = 


P a2L  „ 
6K 


/f(Y)  -2- 

0 

7 f(x)  80 

/02+l  )(3/32+l ) 

_ojl 


302  + 1 


(3.29) 


with  the  retarded  angle  /3  defined  by 
03  + 0 = K(D  + b - X). 


(3.  30) 


The  nondimens ional  quantity  b in  (3.  30)  is  the  distance  of  the  blade  span 
axis  from  the  blade  leading  edge.  This  distance  is  usually  taken  as 
b=  1/4. 

It  may  be  worthwhile  to  review  what  has  been  accomplished  by  the 
process  of  integration  by  parts,  the  large  aspect  ratio  (small  c)  approxi- 
mation, and  transonic  scaling.  The  main  points  are: 

(1 ) The  integral  formula  (3.  29)  for  p is  only  a single  integral. 
Double  integration  over  the  entire  planform  has  been  avoided. . 

(2)  The  slope  distribution  f(x)  with  its  singular  behavior  at  a 
round  leading  edge  has  been  eliminated  in  favor  of  the  thickness  distri- 
bution f (x)  which  is  finite  everywhere. 
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(3)  The  time  derivative  in  front  of  the  integral  that  represents  p 
is  taken  inside  the  integral  without  introducing  new  singular  behavior 
in  the  integrand. 

(4)  The  function  of  retarded  angle  j3  that  appears  in  the  integrand 
of  (3.  29)  is  a regular  function  of  — :t  is  not  singular  at  any  chordwise 
station  * . 

(5)  An  acoustic  similarity  parameter  K appears.  The  pressure  p 
depends  on  ambient  density,  sound  speed,  and  thickness  ratio  in  a trivial 
way.  The  dependence  of  p on  aspect  ratio  (or  its  inverse  e)  and  tip 

Mach  number  M when  M is  near  one  is  entirely  concentrated  in  K. 

2 

This  dependence  is  not  so  simple:  p depends  explicitly  on  K in  (3.  29), 
and  there  is  an  implicit  dependence  on  K through  the  dependence  of  the 
scaled  retarded  angle  3 on  K. 

(6)  The  form  of  an  acoustic  wave  at  a fixed,  large  radial  position 
R is  a function  of  the  characteristic  variable  ft  defined  by  (3.  22).  This 
variable  is  also  a characteristic  variable  of  the  linear  partial  differen- 
tial equation  that  determines  p.  It  was  therefore  clear  from  the  begin- 
ning that  the  quantity  Rp  should  depend  on  R,  6,  and  t through  a variable 
like  ft.  The  form  of  this  variable  makes  explicit  the  scale  width  of  the 
wave  profile.  For  example,  if  the  observer  is  in  the  plane  of  rotation, 
then  at  a given  (large)  distance  r and  angular  coordinate  O' , the  tem- 
poral duration  of  the  wave  — the  time  interval  during  which  there  is  a 
significant,  audible  and  measurable  wave  amplitude  — is  of  the  order  of 
the  parameter  CM;  or,  for  M near  one,  the  width  is  the  order  e,  the 
inverse  of  the  blade  aspect  ratio. 

One  must  not  place  too  much  emphasis  on  any  of  the  items  listed 
above.  Taken  together,  however,  they  produce  a remarkably  simple 
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formula  for  acoustic  pressure.  One  advantage  of  the  formula  is  that  it 
is  easy  to  use  in  a numerical  calculation. 

The  formula  (3.  29)  for  p may  be  written  in  compact  form  as 


» • *A  k F‘n»- 


(3. 31) 


where  the  nondimens ional  profile  of  an  acoustic  wave  is 


F(S1) 


•-4  /«*>*'• 

o L 


<32+1)<3J32+1)J  3/32  + 1 


(3.32) 


The  form  (3.  31  ) for  p is  especially  convenient  when  nonlinear  behavior 
and  acoustic  shocks  are  considered.  One  reason  for  getting  the  form 
(3.  31),  with  the  relatively  simple  form  for  F(f2),  is  a reduction  in  the 
amount  of  numerical  work  required  in  a nonlinearization  scheme.  The 
function  F would  be  identified  in  sonic  boom  work  as  a Whitham 


F -function. 


-L.  fll 
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4.  ACOUSTIC  SHOCKS 


Acoustic  pressure  in  the  linear  theory  appears  in  Section  3 as 

p=  p0^tF(S!> 

where  the  characteristic  variable  £2  is 

wt-01  - wR/a 

£2  = — ° . 

€ M cos 


(4.1). 


(4.2) 


The  profile  of  the  wave  is  F (£2 ), defined  by  equation  (3.  32).  This  profile 
remains  undistorted  with  change  in  £2  , and  the  linear  p is  always  shock 
free.  It  is  known  from  experimental  work  (ref.  8 ) that  sound  waves 
generated  by  high  speed  rotors  show  a sharp,  shock-like  discontinuous 
form  in  the  plane  of  rotation.  These  discontinuous  wave  forms  occur  even 
when  the  blade  is  operated  at  near  zero  lift,  which  suggests  that  the  signal 
and  its  shape  are  caused  by  blade  thickness  and  a slowly  varying  sound 
speed  in  the  far  field.  Inspection  of  wave  forms  calculated  from  the 
simplest  linear  theory  (as  in  Figure  4 ) immediately  show  that  these 
forms  will  roll  over  and  form  a shock,  and  the  shock  may  be  calculated  by 
correcting  the  far  field  sound  speed.  Whether  the  shock  will  have  signifi- 
cant strenght  (strong  enough,  for  example,  to  rattle  windows  at  some 
distance  from  the  blade  system)  is  not  immediately  evident  from  a graph 
such  as  Figure  4 . 

Departure  of  the  local  sound  speed  from  its  ambient  value  aQ  may  be 
found  from  the  nonlinear  partial  differential  equation  for  the  velocity 
potential.  The  asymptotic  form  of  that  equation  at  large  distance  from  the 
blade  surface  is  equation  (2.12) 
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W *99  ' Zu*R^R9  ~ 


ao+  (\-D 


(^rr  + r*  *r) 


(4.3) 


Nonlinearization  of  the  linear  p of  (4.  1)  proceeds  by  finding  a correction 
to  the  constant  sound  speed  aQ.  This  correction  leads  to  improved 
characteristics  for  equation  (4.3)  that  will  replace  the  linear  character- 
istic variable  . 

Let  a characteristic  surface  of  equation  (4.  3)  be  denoted  by 
S(R,  9)  = 0 , and  write  S as 


S = 9 + g(R)  . (4.  4) 

The  surface  S = 0 is  a solution  to  the  eikonal  equation  corresponding 
to  (4.3), 

“2se-2“^sRse  = [ao  + lv-i)«6]s2  . 

Replacing  S by  its  form  (4.  4)  gives  an  equation  for  g (R): 


go 


2u<f>Rg(R)=  [ ao  + (V-D«00j  8 (R). 


(4.  5) 


Treat  0^  and  0^  in  (4.  5)  as  small  and  solve  the  quadratic  equation  for 
g to  first  order  in  these  derivatives  to  get 


(4.6) 


and  0g  in  (4.  6)  by  their 
linear  forms  derived  from  equation  (4.1)  . Pressure  p and  potential  0 
are  related  by 

P = - Po  ^ ' 


The  nonlinearization  scheme  replaces  0^ 
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and  so,  using  (4.  1)  , 

La2 

\ 1—  F (SI)  . 


(4.7) 


The  derivative  is  also 

*^£7  € M cos  4*  ^ £7  * 

where  (4.2)  has  been  used  for  £7^  . Comparing  (4.7) 
the  important  derivative 


x € M C08  4* 

Z 


La' 


F (£7)  . 


(4.8) 

and  (4.  8)  gives 


(4.9) 


The  derivatives  ^ and  <pQ  that  are  needed  in  (4.6)  are  calculated 
from  (4.  9),  retaining  only  radiative  terms  that  are  0(R‘  *)  . We  have 


= 4*7  UR  • $0  = ^£7  • 

and  therefore 

**.  * -V  F <D>  • \ * V1  F <°>  • 

Substituting  these  relations  into  (4.  6)  for  g gives 


g (R)  = 


mL  v+  1 F(£7) 

a0  2 R 


(4.10) 


The  essential  approximation  in  nonlinearization  is  to  treat  £7  in  equa- 
tion (4.  10)  as  constant  and  then  simply  to  integrate  (4.  10)  with  respect 
to  R.  Doing  so  gives 


g(R>  = 


fcp  R 

a. 


X±J_ 

2 


M F (£7)  log 


fe) 


+ c (£7)  , 


(4.11) 
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where  c (ft)  is  the  constant  of  integration  which,  as  indicated,  may 
depend  on  ft  . The  constant  Rq  in  equation  (4.  11)  is  some  radial 
position  at  which  g(R)  would  take  on  its  linear  form  g(R)  =u»R/a^  + c(ft). 
The  assumption  is  that  there  is  some  region  relatively  near  the  radiating 
body  within  which  the  acoustics  may  be  treated  in  a linear  fashion;  this 
region  is  sufficiently  close  to  the  body  so  that  waves  have  not  propagated 
far  enough  to  break,  roll  over,  and  form  a shock  wave.  In  the  helicopter 
application,  if  the  blade  tip  Mach  number  is  sufficiently  small,  there 
will  certainly  exist  a region  of  linear  acoustics,  and  distortion  in  the 
profile  will  truly  be  limited  to  the  far  field.  However,  small  tip  Mach 
number  generally  implies  small  amplitude  acoustics,  and  an  acoustic 
shock  may  be  so  weak  as  not  to  merit  the  name  shock  wave.  On  the  other 
hand,  when  tip  Mach  number  is  near  one  and  the  amplitude  of  an  acoustic 
disturbance  is  relatively  large,  the  acoustic  wave  may  merge  near  the 
blade  with  an  aerodynamic  shock  and  the  whole  idea  of  some  initial,  near 
field,  linear  region  that  supports  the  far  field  acoustics  may  be  lost.  The 
experimental  work  of  reference  8 suggests  that  the  latter  state  of 
affairs  may  indeed  be  typical  of  a high  speed  rotor.  The  authors  of 
reference  8 are  continuing  an  experimental  program  in  which  the  pressure 
field  of  a rotor  is  mapped  in  a continuous  fashion  from  the  blade  surface 
to  the  far  field. 

The  nonlinear  program  described  here  has  deficiencies  (which  will 
be  pointed  out)  but  it  should  be  kept  in  mind  that  a complete,  thorough- 
going nonlinear  calculation  - near  and  far  field  - will  certainly  have 
embedded  in  its  solution  a nonlinear  acoustics  of  the  type  that  we  are 
outlining.  The  crucial  problem  in  a more  accurate  solution  is  to  find  the 
correct  F (ft)  on  which  to  base  the  far  field  solution. 
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Returning  to  equations  (4.  4)  and  (4.  11),  the  characeristic  surface  S 

is 

S = e + ^ MF(0)  log 

o \ *o 

Recall  that  0 is  the  cylindrical  angle  in  coordinates  fixed  to  the  blade.  It 

is  related  to  the  inertial  reference  frame  O'  and  t by  0 = O'-  wt  . Setting 

6 ~ 0'-ut  in  (4.  12)  and  rearranging  terms  gives  an  expression  for  the 

improved  12  - characteristics 


) 


+ c (12)  = 0. 


(4.  12) 


«t-e«-  MF(12)  log(  ) + c(Q)  . (4.  13) 

o \ o / 

The  correctly  scaled  characteristic  variable  12  in  the  linear  theory  of 
Section  3 is 

ut-0'-uR/a 

w --1°  = 12  . 

€ Mcos4< 

Dividing  (4.  13)  by  fM  cos  + and  making  an  obvious  choice  for  c(12)  gives 


wt-fl'-wR/a 


e M cos 


111 
2 e cos 


F(12)  log  . 


(4.  14) 


The  constant  Rq  in  the  argument  of  the  logarithm  in  equation  (4.  1 3)  has 
been  replaced  by  the  blade  span  L . The  choice  of  a length  scale  with 
which  to  nondimens ionalize  the  radial  coordinate  R in  (4.  13)  is  not  obvious. 
Lengths  that  are  or  are  not  suitable  are  based  on  the  following  observations: 

1.  The  acoustic  pressure  is  an  outer  solution  to  the  flow  about  the 
blade.  The  blade  surface  has  three  characteristic  and  disparate  length 
scales:  the  blade  span  L,  blade  chord  I , and  section  thickness  t , 
with  L>>|>>t.  It  is  a general  result  that  the  smallest  length  scales 
do  not  play  a role  in  the  outer  solution  (the  smallest  lengths  are,  so  to 
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speak,  lost  from  our  distant  view  of  the  blade).  The  only  character- 
istic geometric  blade  length  to  survive  in  the  outer  solution  would  be 
the  span,  so  it  might  be  used  in  (4.  13). 

2.  A wave  length  based  on  the  blade  fundamental  frequency 
f = w/27r  is  not  a suitable  scale.  The  entire  acoustic  wave  form  will 
propagate  past  a fixed  observer  with  repetition  frequency  f , but  this 
frequency  is  not  a typical  frequency  in  the  wave  form  itself.  Referring 
to  the  definition  of  the  variable  12  in  equation  (4.  2)  it  is  apparent  that  a 
reference  frequency  for  the  wave  form  is  w/c  M,  or,  for  M near  one, 
the  blade  fundamental  frequency  times  the  blade  aspect  ratio,  which  is  of 
course  much  larger  than  <*»  . Most  of  the  acoustic  energy  is  lodged  in 
frequencies  near  «/e  M.  These  analytical  remarks  agree,  incidentally, 
with  experimental  work;  a significant  part  of  high  speed  helicopter  sound 
is  high  frequency  sound.  Now  let  X be  a wave  length  based  on  the  blade 
fundamental  frequency  w , and  let  X'  be  a wave  length  based  on  the 
angular  frequency  w 1 = w/c  M.  Then 

X'  = € M X . 

Noting  that 


we  have  for  a length  scale  in  the  far  field 
X ' = 2 7T  e L . 

Although  6 is  small,  as  a practical  matter  the  quantity  2 ir  e is  not  so 

small.  For  an  aspect  ratio  of  13.71  (used  in  ref.  ( 8 ) and  in  our 

/ 

calculations)  we  have  X1  = 0.  46L,  and  the  distinction  between  the  far 
field  scales  X'  and  L is  not  a serious  one.  The  simplest  choice  here 


is  L . 
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Now  define  a new  variable  T as 
wt-0'-wR/a 

T = — . (4.15) 

eMcos  ^ 

Equation  (4.  14)  for  the  new  characteristic  variable  ft  is 

T=n-  2 <£'0.4.  f(n)l„6(£.)  . (4.. 6) 


Acoustic  pressure  p in  (4.  1 ) is  a function  of  ft  in  linear  theory  and  is 
the  same  function  of  ft  in  the  nonlinear  theory.  What  is  required  is  p as 
a function  of  T . It  is  therefore  necessary  to  invert  (4.  16)  to  find 
ft  sft  (T,R)  and  substitute  this  result  into  F in  (4.  1).  Equation  (4.  16) 
is  transcendental,  and  so  the  ft  (T,R)  inversion  must  be  a numerical 
calculation.  We  then  have  two  sets  of  results: 

Linear  theory  p = p a*"  — F (ft)  , (4.  17a) 

O O XV 

T = ft  . (4.  17b) 

Nonlinear  theory  p = pQa*  F (ft)  , (4.  18*) 

T * 0 - irhn  FW  l0«  (£•)  • «4  l8b> 

A shock  will  form  when  the  T(0)  relation  (4.  18b)  forms  an  envelope  with 
respect  to  ft  . This  envelope  is  found  by  differentiating  (4.  18b)  with 
respect  to  0 to  give 

1 - 27mVo.->  f <n>  M r)  • 

The  radial  position  R/L  where  the  envelope  first  forms  and  the  shock 

t 

begins  is  determined  by  the  value  of  ft  that  makes  F a maximum.  Thus 


108  (xv 


MIN 


2 € COB  ^ 

y+1 


{ F (ft)] 


MAX 


(4.19) 
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A difficulty  appears:  the  value  of  R/L  predicted  by  (4.  19)  for  a typical 
rotor  and  flight  conditions  is  in  poor  agreement  with  experiment.  Let  us 
use  numbers  in  (4.  19)  for  which  experimental  data  are  available  (ref.  8)  . 
Take  as  aspect  ratio  13.71,  so  that  c = 0.07  3.  Set  4*  = 0 for  an  in-plane 
observer,  and  use  y = 1 . 4 for  air.  Let  the  blade  have  rectangular  plan- 
form  with  NACA-0012  sections.  For  a tip  Mach  number  0.92,  the  theory 
of  Section  3 gives  0.0179  as  the  maximum  value  of  F . Equation  (4.  19) 
then  gives  R/L  = 30  for  the  in-plane  radius  at  which  a shock  begins. 
According  to  the  work  of  reference  8 , at  these  conditions  a shock 
begins  on  the  blade  surface  (the  usual  transonic  aerodynamic  shock,  much 
as  for  a fixed  wing),  but  it  then  continues  off  the  blade  tip  and  into  the  far 
field,  without  interruption.  This  shock  configuration  will  appear  at  tip 
Mach  numbers  less  than  0.  92  for  the  specified  blade  geometry.  This 
discrepancy  between  theory  and  experiment  is  probably  due  to  nonlinear 
effects  that  occur  on  the  blade  surface  near  the  tip,  and  which  persist  for 
some  few  chords  off  the  tip  where  more  distinctively  acoustic  effects  will 
begin.  The  blade  surface  nonlinearities  cannot,  of  course,  be  accounted 
for  by  a simple  modification  of  linear  theory. 

The  above  mentioned  inaccuracy  in  the  nonlinearization  may  be 
partially  remedied.  The  linear  result  has  an  acoustic  pressure  and  an 
F-function  that  are  proportional  to  blade  thickness  6 . It  is  known,  however, 
that  near  the  blade  surface,  the  pressure  distribution  that  is  based  on  the 
nonlinear  differential  equation  (2.  11)  in  proportional  to  6^^  , not  6 . 

This  is  a transonic  scale  effect  that  appears  in  rotor  theory,  just  as  it  does 

in  fixed  wing  theory.  We  may  then  modify  the  F-function  in  equations  (4.  18) 

l/3  2/^ 

by  dividing  it  by  6 , thereby  making  F proportional  to  6 ' . The  result 


is  a 
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Modified  nonlinear  theory; 


2 L F(S2) 

P=  Poao  X 7l7T 

T=«-  ^log  (£ 

2c  cos  4*  61/3 


(4.  20a) 
(4.  20b) 


Equation  (4.  18b)  for  the  shock  initial  radius  is  changed  to 


MIN 


2 € cos  4^ 
V+i 


1/3 


tF  <n>Wx 


(4.21) 


This  rescaling  of  the  F-function  is  not  entirely  arbitrary.  An  exact 
solution  to  the  differential  equation  (2.  3)  would  yield  the  entire  flow  field, 
both  aerodynamic  and  acoustic.  In  the  far  field,  an  approximation  to  that 
solution  would  take  the  form  (4.  20),  and  the  approximation  would  be  a 
good  one.  The  F-function  that  approximates  the  exact  solution  would  not, 
however,  be  the  linear  F used  in  (4.20);  it  would  be  some  other  F that 
approximates  the  far  field  behavior  to  the  numerical  solution  to  the  dif- 
ferential equation  . (We  leave  aside  here  the  fact  that  if  an  exact  solution 
to  the  differential  equation  were  known,  then  a boom  approximation  to  that 
solution  would  not  be  necessary  ; but  this  is  not  to  say  that  such  an  ap- 
proximation would  be  of  no  interest).  The  numerically  derived  F would 
have  the  scaling  that  goes  with  the  differential  equation.  Nevertheless,  it 
has  to  be  acknowledged  that  merely  rescaling  the  nonlinearized  formula 
(4.  18)  to  fit  a known  property  of  an  exact  solution  is  something  of  an  ad  hoc 
adjustment. 

Let  us  find  the  initial,  in- plane  location  of  the  shock  from  the  modified 
formula  (4.  20),  using  the  same  data  that  were  used  in  the  unmodified  for- 
mula (4.18)  ; aspect  ratio  = 13.71,  M « 0.92,  rectangular  planform, 
NACA-0012  sections.  Then  (4.21)  gives  R/L  = 5,  as  opposed 


previous  R/L  = 30.  Movement  of  the  shock  toward  the  tip  is  a desirable 
trend  in  view  of  the  results  of  reference  8.  The  boom  formulas  should,  of 
course,  be  applied  in  the  region  R/L  >>  1 . Their  application  here  is 
treading  into  the  region  R/L  = 5,  which  scarcely  qualifies  as  very  large. 

The  remaining  acoustic  shock  calculations  of  this  report  are  based 
on  the  modified  scheme  (4.  20).  A discussion  of  numerical  results  based 
on  this  and  the  previous  section  is  taken  up  in  Secion  5. 
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5.  NUMERICAL  RESULTS  AND  CONCLUSIONS 

Numerical  results  shown  in  Figures  4-10  are  presented  as  plots 
of  acoustic  pressure  against  a characteristic  variable  at  some  fixed  in- 
plane radius.  Since  the  radius  is  fixed,  the  characteristic  variable 
may  be  considered  a scaled  time  variable.  Linear  calculations  are 
based  on  equation  4.17  ; the  characteristic  dependent  variable  is  S2 
defined  by  equation  (4.  2).  Nonlinear  p is  calculated  using  (4.  20a)  with  a 
characteristic  dependent  variable  T defined  by  (4.20b).  In  the  linear 
theory,  and  T are  identical;  in  a nonlinear  calculation,  O is 
eliminated  in  favor  of  T , the  functional  relation  being  equation  (4.  20b). 

Some  of  the  blade  geometry  is  fixed  in  Figures  4-10:  planform 
is  rectangular,  there  is  no  twist  or  taper  in  thickness,  thickness  ratio  is 
12%  , aspect  ratio  is  13.71.  All  calculations  are  for  hover.  Tip  Mach 
number  is  0.  92  and  pressure  is  calculated  in  the  plane  of  rotation.  Two 
blade  sections  are  used,  biconvex  and  NACA-0012.  Geometry  and 
flight  conditions  were  selected  for  comparison  with  experimental  data  of 
reference  8. 

It  is  not  necessary  to  recapitulate  here  the  well  known  results  of 
elementary  sonic  boom  theory  on  which  our  nonlinear  results  are  based. 
The  method  and  elaborations  on  it  appear  in  Whitham's  book.  It  is  enough 
to  note  that  the  system (4.  18)  or  ( 4.  20  ) lead  to  wave  form  distortion  and 
finally  to  wave  breaking  and  shock  formation.  Location  of  a shock,  when 
one  occurs,  and  calculation  of  portions  of  the  wave  profile  separated  by 

‘ k 

the  shock  wave  are  found  by  the  usual  methods. 

Two  features  that  are  currently  used  in  an  initial  comparison  of 
theoretical  results  with  experimental  data  are  the  peak  negative  pressure 
and  the  general  shape  of  the  pressure  distribution  as  a function  of  time. 
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or  as  a function  of  some  equivalent  time  variable  (here  our  12  or  T ). 

A completely  linear  calculation  for  a blade  with  biconvex  sections 
is  shown  in  Figure  4.  The  graph  is  symmetrical  about  12  = 0 because  the 
span  axis  of  the  blade  is  taken  to  be  the  section  half-chord  point  for  bi- 
convex profiles;  the  span  axis  is  therefore  an  axis  of  symmetry  of  the 
blade.  Figure  5 shows  a linear  calculation  for  a blade  with  NACA-0012 
sections.  Flight  conditions  for  Figures  4 and  5 are  the  same.  Figure  5 
also  has  a symmetrical  appearance,  but  here  the  axis  of  symmetry  goes 
through  a value  of  12  slightly  greater  than  zero.  This  is  because  the  span 
axis  of  the  0012  blade  was  taken  as  the  section  quarter-chord  point,  as  it 
would  be  for  a blade  in  actual  flight  conditions. 

A surprising  result  is  that  apart  from  the  location  of  axes  of  sym- 
metry, Figures  4 and  5 are  nearly  identical.  The  peak  negative  pressure 

2 

in  both  cases  is  -481N/m  . According  to  reference  5,  the  measured  peak 
negative  pressure  in  approximately  -960N/m^  for  a l/7  scale  model  of 
the  UH-1H  helicopter  blade  at  the  flight  conditions  of  Figure  5.  The 
measured  peak  is  therefore  larger  in  magnitude  than  the  linear  result  by 
very  nearly  a factor  of  2.  This  discrepancy  seems  to  be  typical  of  all 
linear  calculations. 

A characteristic  of  Figures  4 and  5 that  affects  harmonic  content  is 
the  temporal  scale  width  of  the  pressure  profile.  We  define  this  width  to 
be  the  time  interval  that  separates  the  two  positive  pressure  peaks.  The 
real  time  width  At  is  related  to  the  separation  Aft  of  the  peaks  in 
Figures  4 and  5 by  At  = ( e M/u  ) Aft  (this  follows  from  the  definition 
of  ft  in  equation  (4.  2)  ) . The  angular  frequency  u>  may  be  determined 
from  u>  =(aQ/L.)M.  For  ambient  sound  speed  a 341m/s  and  a blade 
radius  L = 1.  04  m (the  scale  model  radius  of  ref.  5)  ,we  have 
w = 302  rad/s  at  a tip  Mach  number  of  0.  92.  The  separation  Aft  in 
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Figures  4 and  5 is  nearly  2.  The  real  time  width  is  therefore  At  = 0.  4 

milliseconds.  The  corresponding  measured  widths  taken  from  reference  5 

are  0.  36  ms  for  M =0.9,  and  0.  43  ms  for  M = 0.  96.  Our  calculated 

value  of  0.  4 ms  therefore  agrees  quite  well  with  experiment. 

Figures  6 and  7 show  nonlinearized  calculations  for  the  biconvex 

and  0012  blades.  There  is  one  startling  result  in  these  figures:  the  peak 

2 

negative  pressure  is  -975  N/M  , and  this  number  agrees  almost  precisely 

2 

with  the  measured  value  of  about  -960  N/m  in  reference  5.  This  result 
must  be  accepted  with  some  caution  when  we  recall  the  particular  kind  of 
nonlinearization  used  to  generate  Figures  6 and  7.  The  usual,  straight- 
forward sonic  boom  nonlinearization  has  been  modified,  as  already 
explained  in  Section  4,  by  rescaling  the  F-function  to  make  it  proportional 
to  thickness  ratio  to  the  two-thirds  power.  This  rescaling  affects  amplitude, 
location  of  a shock  wave  if  one  occurs,  and  the  amount  of  distortion  in  the 
wave  profile  - it  permeates  the  whole  calculation.  The  modification  is 

suggested  by  observing  that  pressure  in  the  exact,,  nonlinear  transonic 

2/3 

theory  is  proportional  to  6 ' near  the  blade.  The  observation  to  make 
here  is  that  agreement  of  calculated  results  in  Figures  6 and  7 with  ex- 
perimental data  may  or  may  not  be  fortuitous.  Calculations  at  other  hover 
tip  Mach  numbers  are  not  yet  available  for  comparison  with  the  full  range 
of  tip  Mach  number  data  of  reference  5. 

It  will  be  noticed  that  Figures  6 and  7 are  distorted  versions  of 
4 and  5.  The  right-hand  portions  of  6 and  7 are  steeper  than  the  correspond- 
ing parts  in  4 and  5.  In  fact,  the  right-hand  portions  of  6 and  7 are  so  steep 
as  to  suggest  the  appearance  of  a shock  wave,  but  this  is  not  the  case  for 
the  data  shown  in  the  lpgends  of  the  two  figures. 
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Figure  8 is  a graph  of  acoustic  pressure  for  a biconvex  blade,  now 
calculated  at  the  breaking  radius  R/L  = 3.  25  . This  radius  is  determined 
by  equation  (4.  21).  The  steepest  portions  of  Figures  8 and  9 are  very 
steep  indeed:  casual  inspection  of  these  graphs  might  suggest  that  shocks 
occur  here.  Of  course,  a shock  must  appear  in  the  pressure  profile  at 
any  radius  greater  than  the  breaking  radii  of  Figures  8 and  9. 

Figure  10  is  a wave  profile  for  the  biconvex  blade,  calculated  at 
R/L  = 6.  48.  A full-fledged  shock  wave  now  appears.  The  shock  extends 
almost  from  the  peak  negative  pressure  to  the  peak  positive  pressure; 
the  steepest  part  of  the  profile  is  almost  all  shock  wave.  The  physical 
dimensions  of  the  anechoic  facility  used  to  generate  the  data  of  refer- 
ence 5 do  not  allow  measurement  at  radial  distances  corresponding  to 
Figure  10,  when  a 1 meter  radius  blade  is  used,  so  comparison  with 
experiment  cannot  be  made  here. 

The  main  conclusion  to  be  drawn  from  Figures  4-10  is  that  a non- 
linear method  gives  the  distortion,  and  finally  the  shock  wave,  that  brings 
calculated  results  closer  to  experimental  data  than  does  any  linear  theory. 
These  results  are  a forceful  reminder  that  good  nonlinear  technique, 
whether  in  the  Ffowcs  Williams -Hawkings  equation  or  in  the  complete 
rotor  differential  equation,  is  essential  if  the  full  range  of  experimental 
data  is  to  be  accounted  for. 

Several  qualitative  conclusions  emerge  from  a discussion  of  the 
nonlinear  potential  equation  that  governs  the  rotor  near  field  aero- 
dynamics and  the  far  acoustic  field.  Two  nonlinearities  appear.  The 
first  one  is  effective  near  the  blade  surface  and  controls  the  local 
aerodynamics,  including  the  shock  wave  that  may  occur  on  the  blade; 


it  is  not  effective  in  the  far  field.  The  second  nonlinearity  is  decisive 
in  maintaining  distortion  and  a shock  wave  in  the  propagating  acoustic 
profile;  it  has  no  counterpart  in  fixed  wing  (M  < 1)  aerodynamics.  Scal- 
ing arguments  show  this  second  nonlinearity  to  be  unimportant  on  the 
blade  surface. 

The  characteristic  surfaces  of  the  rotor  differential  equation  may 
shed  some  light  on  recent  experimental  results.  It  has  been  observed 
that  a hovering  transonic  blade  has  a kind  of  lower  critical  Mach  num- 
ber at  which  abrupt  changes  in  the  acoustic  field  appear.  At  and  above 
this  Mach  number,  the  blade  surface  and  the  far  field  seem  to  become 
linked  by  an  uninterrupted  shock  wave,  and  the  sound  field  becomes 
measurably  more  intense.  The  appropriate  conjecture  is  that  this 
behavior  is  initiated  when  the  extreme  blade  tip  first  touches  or  punc- 
tures the  sonic  surface  defined  by  the  nonlinear  potential  equation. 

The  validity  and  usefulness  of  this  idea  will  require  numerical  solution, 
of  the  differential  equation  and  more  detailed  experimental  investigation. 
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6.  NOTATION 


a Constant  adiabatic  sound  speed  = 341  m/s  . 

o 

b Distance  of  blade  span  axis  from  blade  leading  edge  as  a fraction  of 

blade  chord. 

D Doppler  factor  = (1-M^)  * . 

f Blade  section  thickness  distribution;  blade  fundamental  frequency. 

F(D)  Profile  of  acoustic  wave. 

i,  j,  k Unit  vectors  along  x,y,  z,  or  ^ > 1 ^ blade  fixed  axes, 

ij  j,'  k'  Unit  vectors  along  x‘,  y',  z'  nonrotating  axes. 

K Acoustic  similarity  parameter. 

I Blade  chord 


Blade  radius  (span). 
Blade  tip  Mach  number. 
Radiation  Mach  number. 


Disturbance  pressure;  acoustic  pressure. 

Cylindrical  radial  distance;  in- plane  radial  coordinate. 

Far  field  unit  radiation  vector. 

Spherical  radial  distance. 

I*-*.!- 

Observer  position  vector. 

Position  vector  of  point  on  blade  surface. 

Unit  radiation  vector. 

S * O : Characteristic  surface  of  rotor  differential  equation. 
Projection  of  blade  surface  onto  plane  of  rotation. 

Time. 

T ime  ( * t ) . 


-LQ- 


T Characteristic  variable  in  nonlinear  theory. 

v Projection  of  blade  velocity  along  outward  normal, 

n 

V Velocity  of  air  relative  to  blade, 

x,  y,  z Blade-fixed  axes, 
x'.x’.z'  Nonrotating  axes. 

o Observer  spherical  angle  measured  from  axis  of  rotation. 

/3  Scaled  retarded  angle, 

y Ratio  of  specific  heats. 

6 Thickness  ratio. 

e Inverse  of  blade  aspect  ratio  I /L  . 

S’  ^2’  S Nondimensional  blade  fixed  coordinates. 

H i»r^2»T^3  Dimensional  blade  fixed  coordinates. 

0 Cylindrical  polar  angle  relative  to  span  axis. 

0'  Cylindrical  polar  angle  in  nonrotating  frame. 

£,  Retarded  angle  = n /2  - (9'-  <*>t* ). 

3 

p Ambient  density  = 1.23  kg/m  . 

o 

T Retarded  time  . 

4>  Velocity  potential. 

4»  Observer  meridian  angle  (out-of-plane  angle)  . 

w Blade  angular  velocity. 

oT  Blade  angular  velocity  vector. 

U Scaled  characteristic  variable. 


Id  inters 


FIGURE  4 
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FIGURE  5 
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FIGURE  6 


Slope  finite  everywhere 


P N/m 


Infinite  slope  at  T 
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